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' Abstract. According to the standard picture for binary interactions, the outcome of binaries surviving the 

evolution through a common envelope (CE) and spiral-in phase is determined by the internal structure of the 
donor star at the onset of the mass transfer, as well as the poorly-known efficiency parameter, t/ce , for the ejection 
of the H-envelope of the donor. In this Research Note we discuss the bifurcation point which separates the ejected, 
unprocessed H-rich material from the inner core region of the donor (the central part of the star which will later 
contract to form a compact object). We demonstrate that the exact location of this point is very important for 
CIh' evaluating the binding energy parameter, A, which is used to determine the post-CE orbital separation. 

Here we compare various methods to define the bifurcation point (core/envelope boundary) of evolved stars with 
masses 4, 7, 10 and 20 Mq. We consider the specific nuclear energy production rate profile, the change in the mass- 
density gradient (Bisscheroux 1998), the inner region containing less than 10% hydrogen, the method suggested 
j^j , by Han et al. (1994) and the entropy profile. We also calculated effective polytropic index profiles. 

The entropy profile method measures the convective boundary (at the onset of flatness in the specific entropy) 
which is not equivalent to the core boundary for RGB stars. Hence, this method is not applicable for RGB stars, 
unless the actual bifurcation point of a CE is located at the bottom of the outer convection zone (resulting in 
larger values of A and larger post-CE orbital separations). On the AGB, where highly degenerate and condensed 
cores are formed, we find good agreement between the various methods, except for massive (~ 20 Mq) stars. 
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1. Introduction 

The simple formalism introduced by Webbink (1984) and 
de Kool (1990) to estimate the orbital decay of binaries 
evolving through a CE and spiral-in evolution requires 
knowledge about the core mass, M core , inside the bifur- 
cation point of the donor star as well as the A-parameter, 
which is a numerical factor introduced in order to correct 
their simple formula for estimating the binding energy of 
the stellar envelope. It has been common practice to use 
a constant value of A = 0.5 for the A-parameter. We have 
demonstrated (Dewi & Tauris 2000) that this is not a good 
approach, since A varies strongly throughout the evolution 
of a star. We concluded that all observations of binary pul- 
sars originating from a CE evolution are consistent with 
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Vcf, < 1, as a result of large values of A possible when the 
internal thermodynamic energy is included. 

We also raised the important question of how to define 
the core mass boundary - i.e. once the spiral-in process 
is initiated, and frictional torques deposit kinetic energy 
in the envelope, where is the exact location of the point 
of bifurcation in the envelope which separates the ejected 
material from the remaining condensed core region ? This 
question is the key issue in this Research Note. It is neces- 
sary to know M core accurately, since it is used to determine 
the exact value of A, which is a strongly increasing func- 
tion exactly at the core/envelope transition. Hence, the 
estimates of the post-CE orbital separation also depend 
strongly on the value of M core . 

We refer the reader to Dewi & Tauris (2000) and ref- 
erences therein for further details on the topic and the 
relevant energy equations and parameters. 
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2. Structure profiles of evolved stars 

In Fig. 1 we show an example of our calculations of the 
stellar structure for a 10 M© star. The results for a 4, 7 
and 20 M Q star are summarized in Table 1. We used a 
chemical composition of X = 0.70 and Z = 0.02, and a 
mixing- length parameter of a = l/H p = 2.0. Convective 
overshooting was also taken into account. Here we used an 
overshooting constant of S ov = 0.10 (Pols et al. 1998). We 
obtained almost the same core masses using S ov = 0.12. 
For the 10 and 20 M© star we assumed a wind mass- 
loss rate according to de Jager (de Jager et al. 1988; 
Nieuwenhuijzen & de Jager 1990). Hence, at the tip of the 
AGB these two stars are estimated to have evaporated 
down to 9.59 and 16.2 M Q , respectively. We assumed no 
mass loss from the 4 and 7M stars. 

We will now briefly discuss the calculations represented 
in Fig. 1. In the bottom panels the two solid lines show 
our calculations of A. We distinguish Ab (calculated from 
the total binding energy) from A g (as derived from gravi- 
tational binding energy alone) - sec Dewi & Tauris (2000). 

2.1. Energy production rate and chemical composition 

A very crude definition of the stellar core region is sim- 
ply as the inner part containing less than a certain mass 
fraction of hydrogen. As an example, we used X < 0.10 in 
Dewi & Tauris (2000) to define the core of the donor star. 
This is shown in the top panels of Fig. 1 as the area to 
the left of the gray shaded region. The dotted line is the 
core radius as a function of its mass. 
A somewhat more physically meaningful approach (al- 
though not necessarily more correct in practice) is to de- 
fine the core region below the (outer) shell burning source 
- for example, below the point of maximum energy pro- 
duction inside the H-burning shell. In the top panels the 
solid line shows the local nuclear energy production rate. 
The shell burning regions are clearly seen. The most abun- 
dant local chemical elements are written in the panel. As 
can be seen from the figure and Table 1, this definition 
sometimes yields significantly lower values of M COIC . 

2.2. The binding energy profile 

The solid lines in the second panel are sinh _1 (AVF), where 
AW = E glav + Ei nt or AW — E glav is the binding energy 
of the envelope to the core (with or without the inclu- 
sion of the internal thermodynamic energy, respectively). 
Han et al. (1994) introduced this function to define the 
core mass boundary. The location is assumed to be at the 
transition between the strongly increasing AW and the 
outer region where AW varies slowly with mass. We used 
the intersection of straight-line fits to find M coro . 

2.3. The mass-density gradient 

The dotted line in the second panel shows the 
mass-density profile. The arrow indicates the point 



where d 2 \ogp/dm 2 = 0. This criterion was used by 
Bisscheroux (1998) to locate the core boundary. The dis- 
advantage of this method, from a technical point of view, 
is that there is not always a unique solution to the equa- 
tion (for small mass steps). 

2.4. The entropy profile and effective polytropic index 

The entropy profile for the 10 M star is plotted as a solid 
line in the third panel. Here "entropy" refers to the lo- 
cal specific entropy per baryon in units of the Boltzmann 
constant. The arrow indicates the sharp onset of the flat 
entropy gradient - another criterion for determining the 
bifurcation point. The core masses at the tip of the RGB, 
defined by this method, are always larger than those de- 
rived from other criteria. The clear discrepancy between 
the entropy profile method and the other alternatives on 
the RGB, leading to different estimates of the remaining 
mass after spiral-in, is easy to understand. The bifurca- 
tion point obtained from the entropy method is expected 
to be located further out, since it is based on the transition 
between the convective and the radiative layer in the stel- 
lar envelope. Unlike evolved low-mass (1M ) stars, these 
more massive stars do not have an outer convection zone 
which penetrates all the way down to the H-shell near 
the standard core boundary region (e.g. Kippenhahn & 
Weigert 1990). Hence, the entropy method results in (too) 
large a core mass on the RGB. 

The dotted line in the third panel is the "effective poly- 
tropic index" defined by: n = 1/(7 — 1) with the "adia- 
batic" index: 7 = d In P/d In p. Giant stars can be consid- 
ered as a condensed polytrope with a core point mass sur- 
rounded by an extended isentropic n — 3/2, 7 = 5/3 enve- 
lope. Less evolved stars are more realistically described as 
composite polytropes consisting of n — 3, 7 = 4/3 cores 
with n = 3/2,7 = 5/3 envelopes (Hjellming & Webbink 
1987). Here, we are mostly interested in the clear discon- 
tinuous behaviour of n at the core boundary in order to de- 
termine its exact location. This discontinuity represents an 
important transition in the stellar structure and could be 
interpreted as the critical point of bifurcation at the end of 
the spiral- in and envelope ejection phase (cf. Section 3.2). 

3. Discussion 

3.1. Mass transfer initiated on the RGB vs. AGB 

We have now demonstrated how different methods result 
in different core masses on the RGB and hence different 
values of A. As an example, one can consider a 4 M star 
at the tip of the RGB. Here the maximum energy produc- 
tion rate and the 10% hydrogen criterion yield core masses 
of 0.58 and 0.59 M , respectively. The mass-density gra- 
dient method results in M coro = 0.68 M and the outcome 
of applying Han et al.'s method is approximately 0.64 M . 
Using the entropy profile and the effective polytropic in- 
dex method yields M corc = 0.74 M . The spread in M core 
from the different approaches results in: 0.6 < Ab < 1.8 
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Fig. 1. The internal structure of a 10 M star at the tip of the RGB (left) and at the tip of the AGB (right). See text. 
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Table 1. The values of M core /M Q , A g and Ab estimated from calculations of four different stars - see text for discussion. 



method 



max e nuc 

X < 0.10 

d 2 log p/ dm 2 — 

Han et al. 

entropy profile 



M = 4.0 M(T 



© 



tip of RGB 
R=67R 



M c . 



0.58 
0.59 
0.68 
0.64 
0.74 



0.32 0.62 

0.36 0.70 

0.66 1.60 

0.59 1.17 

0.73 1.81 



tip of AGB 
R=1040 R 



Mc, 



1.37 
1.37 
1.37 
1.37 
1.37 



0.91 
0.91 
0.91 
0.91 
0.91 



M = 7.0 M, 



o 



AGB 

R=374R 



A E Ab M c , 



A B Ah 



1.58 
1.80 
1.80 
1.80 
1.80 



0.10 0.20 
0.59 2.85 
0.59 2.85 
0.59 2.85 
0.59 2.85 



M = 10.0 Mp 



© 



tip of RGB 
R=374R 



Mc, 



A E Ah 



1.99 
2.08 
2.54 
2.52 
2.62 



0.09 0.13 
0.11 0.21 
0.54 1.80 
0.51 1.75 
0.60 2.23 



tip of AGB 
R=588R 



Mc, 



A s At 



2.47 
2.84 
2.85 
2.85 
2.86 



0.06 0.10 
0.40 1.80 
0.46 2.20 
0.46 2.20 
0.55 3.20 



M = 20.0 M f 



o 



tip of AGB 
R=1040 R 



Mc, 



Ag Ah 



6.80 
6.80 
7.50 
7.60 
7.80 



0.02 0.05 
0.02 0.05 
0.22 0.95 
0.32 1.70 
0.45 3.50 



and 0.3 < A g < 0.75 (with and without internal thermo- 
dynamic energy). If we now assume the simple CE en- 
ergy equations (Webbink 1984; de Kool 1990) to be ap- 
proximately valid, and this 4M© (R = 67 R Q ) star is in 
a common envelope with an in-spiralling 1.3 M© neutron 
star, the evolution would then result in a final post-CE or- 
bital separation of 0.70 or 2.75 R©, for M coro = 0.58 and 
0.74M©, respectively. Here we assumed A = Ab, JJce = 1-0 
and no significant accretion onto the neutron star during 
spiral-in. However, the core radius of the stripped star is 
0.1 R© (M corc = 0.58 Al©) or 3.5 R© (M corc = 0.74 M©). 
Hence, in the latter case, the large core would then have 
a radius larger than the predicted post-CE separation. 
Therefore this system would coalesce during the spiral-in 
- perhaps leaving behind a black hole. Assuming A = A g 
(i.e. considering gravitational binding energy alone), or 
choosing a small efficiency parameter (e.g. rycE = 0.3), 
would decrease the predicted post-CE separation still fur- 
ther and lead to an earlier coalescense. We conclude that 
knowledge of the accurate bifurcation boundary is crucial 
to forecast the outcome of a CE evolution on the RGB. 

Many binaries come into contact while the donor star 
is expanding enormously on the AGB. If we consider the 
same 4 M© star as before, but now at the tip of the AGB, 
all methods yield exactly the same core mass. For the 7 
and 10 M© stars on the AGB, the different methods also 
yield the same core mass (except for the location of the 
maximum energy production rate in the H-shell, which is 
below this boundary at the outer edge of the shell source). 
It is important to notice the fact that the envelope binding 
energy is a strongly increasing quantity near/at the core 
boundary on the AGB. Hence, the resulting values of A 
are also increasing sharply at this transition (see Fig. 1). 
Since the material of the donor star is peeled off from the 
outside, we suggest using the (larger) value of lambda at 
the mesh point, just outside the bifurcation point. 

3.2. Envelope ejection process 

While the location of the bifurcation point, separating 
the ejected envelope from the condensed core region, is 
straightforward to define on the AGB, it is important to 
know for RGB stars if the envelope of the evolved donor is 



ejected at the bottom of the convection zone, or closer to 
the stellar center at the core boundary near/at the outer 
shell source (surrounding the inner part of the star). For 
convective material in an isentropic envelope, R oc M -1 / 3 
and thus removal of envelope mass causes the star to ex- 
pand. Therefore, in a CE system, the envelope is easily 
peeled off until the entropy profile increases outwards (in 
the deeper radiative layers). From this point on the star 
will no longer expand and the envelope may already be 
separated from the inner region, depending on the details 
of the ejection process. It would be interesting to see if 
multi-dimensional hydrodynamical calculations could an- 
swer this question. 

Binaries which come into Roche-lobe contact early in the 
stellar evolution of the donor will not survive the spiral-in 
process at all (there is simply not enough orbital energy 
available to expel the envelope). 

The values of A presented in Table 1 in Dewi & Tauris 
(2000) were calculated from the 10% hydrogen criterion. 
From the work presented here we can see that these val- 
ues should therefore be taken as minimum values of A, 
supporting our previous conclusion that tjce need not be 
large since the value of A itself is expected to be large. 
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